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ABSTRACT 



Motivated by the desire to put a free electron laser (FEL) weapon on a ship, 
the FEL and the related process of Compton backscattering are studied. The theme of 
the majority of this work is the interaction of the Gaussian optical mode with a beam of 
relativistic electrons. 

Classical FEL theory is reviewed in Chapter II. Simulations based on the 
classical theory are used in Chapter III to study a proposed 1 kW (kilowatt) infrared 
FEL. In Chapter IV, simulation is used to study the problem of electron beam/optical 
mode overlap in an ultraviolet (UV) FEL. A new concept, the FEL with a short 
Rayleigh length, is studied in Chapter V. The idea is tested on the UV FEL, then used 
to design and simulate a megawatt-class FEL for ship self-defense. 

An analytical calculation of the Compton backscattering of laser light is 
performed in Chapter VI. A quantum electrodynamics (QED) formalism is used to find 
the spectrum and angular distribution of photons scattered out of a Gaussian optical 
mode by relativistic electrons. 
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I. INTRODUCTION 



Today’s US Navy is the most powerful the world has ever known. It can 
position aircraft, ships, and submarines within striking distance and project its power 
deep within the borders of almost any country. Yet today’s ships are arguably more 
vulnerable to attack from mines and anti-ship cruise missiles (ASCMs) than ever 
before, as a result of the Navy’s shift from blue to brown water. Ship defense against 
ASCMs is mainly a problem of reaction time; that is, there is not much time (usually 
seconds, depending on the missile’s speed and the ship’s detection range) between 
detection of the missile and impact of the missile with the ship. In this case a self- 
defense weapon which travels at the speed of light would be most useful. With 
continuing improvements in ASCM speed, stealth, and maneuverability, the problem 
gets even more difficult for existing defenses. A directed energy weapon (DEW) could 
be the solution. 

Experiments in missile destruction and atmospheric propagation have shown 
that a megawatt (MW) class free electron laser (FEL) is the directed energy weapon of 
choice for ship self defense (SSD). Megawatts are needed in order to deposit the 
required energy on the target in a short amount of time. More power probably would 
not shorten this time because thermal blooming, which results from atmospheric 
absorption, will degrade propagation. The FEL is needed because of its inherent 
tunability and designability, characteristics no other laser system can boast. The laser 
weapon’s wavelength must operate in very narrow atmospheric absorption windows in 
order to propagate the beam to the target, avoiding problems like thermal blooming. 
The FEL has another advantage over other laser systems in that its gain medium is a 
vacuum. 

While SSD is seen as the mission area where laser weapons can get their 
start, a MW-class FEL aboard a ship could be an awesome offensive weapon as well. 
It was demonstrated in the 1970s that it is easier for the laser weapon to destroy 
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crossing targets than incoming targets. It is possible that the precision and power of 
the FEL could also be used against aircraft and their pilots, satellites, ballistic missiles, 
other ships, and even personnel. Unlike the proposed Arsenal Ship, which merely 
packs a ship with more missiles, a shipboard directed energy weapon would truly bring 
about a revolution in military affairs, which would hasten the obsolescence of the 
missile. In this work we study the FEL and related phenomena, contributing to a field 
of scientific research with the goal of putting a DEW on a ship. 

Several topics are presented on the interaction of laser beams with relativistic 
electrons; this includes the free electron laser. In the next chapter, an overview of 
classical FEL theory necessary for understanding the ensuing chapters is provided. 
We describe the components of the FEL and then find the emission wavelength using 
basic physical arguments. Using the Lorentz force and Maxwell’s equations, we derive 
a set of two equations called the FEL pendulum and wave equations which describe 
the FEL interaction. Using these equations, we find expressions and describe the 
mechanisms for FEL gain. Modifications of the FEL interaction due to short optical 
pulses or the Gaussian optical mode are considered in the last two sections of the 
chapter. 

Original work simulating a proposed one kilowatt (IkW) upgrade to the 

FIREFLY FEL at Stanford University is presented in Chapter III. This upgrade would 
represent a one hundred-fold increase in FEL average power over that which currently 
exists. This work was originally presented at the 16th International FEL Conference 
and published as "Simulations of the Stanford FIREFLY IkW free electron laser" in 
Nuclear Instruments and Methods in Physics Research (NIM) A 358, ABS 44. (1995). 

In Chapter IV is presented original research concerning optical mode 
transverse effects in ultraviolet free electron lasers (UVFEL). Specifically, it is the 

UVFEL proposed by the Laser Processing Consortium (LPC). a group based at the 

Thomas Jefferson National Accelerator Facility (TJNAF), that is simulated. The results 
were presented at the 17th International FEL Conference and published as 
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"Transverse effects in UV FELs" in NIM A 375, ABS 61, (1996). 

New and original work on a new concept in FEL design, the short Rayleigh 
length FEL, is presented in Chapter V. The concept is first tested on the LPC UVFEL. 
Those results were presented at the 18th International FEL Conference and will be 
published in 1997. We then use this concept to design and simulate a MW-class FEL 
oscillator suitable for ship self-defense. This design and the simulation results were 
presented at a workshop at TJNAF entitled "Navy MW-Class SSD FEL Concepts" 
which was sponsored by the Navy High Energy Laser office at the Space and Naval 
Warfare Command (SPAWAR). 

In Chapter VI we present new and original calculations of Compton 
backscattering of laser light from electrons. This process has been proposed as a 
means of generating x-rays for many applications, including medical imaging, and y- 
rays for a so-called y-y collider. Our calculations are the first attempt to incorporate the 
finite size of the laser beam into the usual scattering problem. 
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II. CLASSICAL FREE ELECTRON LASER THEORY 



A. THE FREE ELECTRON LASER 

The Free Electron Laser is a device which transforms the kinetic energy of a 
beam of relativistic electrons into coherent radiation (Ref. 1). A transverse, static 
magnetic field causes the electrons to wiggle and thereby radiate electromagnetic 
energy (Ref. 1). The radiation can either be stored in an optical cavity and amplified 
over many passes, as in the FEL oscillator, or amplified in a single pass, as in the FEL 
amplifier (Ref. 1). In this work, we will be concerned only with the oscillator 
configuration. 

The major components of the FEL oscillator are the electron accelerator, 
undulator, and optics. The accelerator delivers a beam of relativistic electrons to the 
FEL with energies anywhere from a few MeV to GeV depending on the wavelength of 
laser light desired. There are many different types of accelerators suitable for use in 
an FEL, but we will only be concerned with the radio-frequency linear accelerator (RF 
LINAC). The RF LINAC accelerates pulses of electrons, called micropulses, using 
either standing or travelling waves of an electromagnetic field (Ref. 2). One 
micropulse can be accelerated within each wavelength of the RF field, called a 
"bucket." For example, if the frequency of the RF radiation is f = 500 MHz, buckets 
are separated by a distance elf =60 cm. Electron micropulses are created and 
inserted into the accelerator by the injector (or electron gun). In our example, if the 
injector also operated at 500 MHz, it would fill every RF bucket and micropulses would 
be separated by 60 cm. The current in the micropulse is called the peak current /. 
and the average current is calculated using the duty factor (the ratio of the micropulse 
length to the distance between buckets). If we take / = 500 A and a micropulse length 
l e = 0.6 mm (l e /c = 2 ps), then the duty factor in our example accelerator is 0.001. 
The average current is then 7 = 500x0.001 A = 0.5 A. 
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The undulator is a set of permanent magnets or electromagnets arranged to 
provide a periodic, transverse, magnetic field of wavelength X 0 along which the 
electrons traverse. The polarization of the undulator’s magnetic field can be either 
helical or linear. The world’s first FEL used a helical undulator, but the linear 
undulator is more common now (Ref. 3). Undulators can range in length from 1m to 
25m. 

The mirrors of the optical cavity (in the oscillator) store the light which is then 
amplified over successive passes of the electrons through the undulator. Total optical 
cavity losses, including the outcoupled light, are quantified by the quality factor Q, so 
that losses per round trip of the optical pulse are (100/0)% (Ref. 3). The optical cavity 
is generally centered on and several times longer than the undulator. 

B. THE RESONANCE CONDITION 

With a basic knowledge of the physics of the FEL interaction we can define 
some important parameters and find the wavelength of the emitted radiation. 

Depicted in Figure 1 is the FEL interaction: the interaction of relativistic 
electrons with the undulator magnetic field, £? u , and the radiated electric, E: r , and 
magnetic, B r , fields. At the top is a schematic of an FEL showing the electron beam 
path through the undulator, the orientation of the undulator magnets, and the mirrors. 
The middle sketch is a blow-up of one wavelength X 0 of the undulator. The undulator 
magnetic field oscillates in the yz-plane, causing the electrons to wiggle in the xz- 
plane. At the bottom of the figure is a blow-up of one half of one wavelength X of the 
radiation field. A portion of the electron trajectory is drawn at an angle to the z-axis. It 
is the transverse component of the electron velocity, that which projects on the xy- 
plane, which couples to the radiation electric field, allowing energy to be transferred 
between them. 
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INCOMING 

ELECTRONS 



OUTGOING 

ELECTRONS 






Figure 1. Relativistic electrons interacting with an 
undulator magnetic field. B u . and radiated electric. B r . 
and magnetic. B r , fields (after Ref. 4). 
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In order for the radiation in the oscillator cavity to grow, the electrons must lose 
energy to the light wave on each pass through the undulator. The Lorentz force 
shows that this energy loss is dy/dt «= $ • where $ = i tic, V is the particle velocity 
and the relativistic factor y= (1 -^ 2 ) -1/2 (Ref. 5). Since is transverse, that is, in a 
plane perpendicular to the z-axis, only a transverse component of couples to the 
cavity radiation. Driven by the undulator magnetic field, that transverse component is 
of the form p x °= R e(Ke lk ° z ) where K is a constant called the "undulator parameter," 
k 0 = 2 n/X 0 , and X 0 is the undulator wavelength. The electric field is of the form 
E r « Re (E 0 e l( ‘ kz ~ a>t+ ^' 1 ) where k - 2n/X, X is the radiation wavelength, to = kc, and <]> is 
an optical phase factor. Therefore, we find that dy/dt « Re (/<E 0 e ,1(/c+ * o)z ~ Q>fl e ,<|) ). It is 
evident from these preliminary considerations that the electron "ponderomotive" phase 
defined as 

C = (k + k 0 )z - (Of (1) 

is fundamental to the FEL interaction. For relativistic FELs, this phase represents the 
position of the electron within a radiation wavelength. 

The time it takes the electrons to traverse the undulator is Lie p z = 10 -8 . 
Rather than work with time scales on the order of 1CT 8 seconds, we introduce a 
dimensionless time x, defined as 

cp 7 ' ct 



\ = 



(2) 



L L ’ 

where L is the undulator length and cP 7 = v z = c in the relativistic limit. As the 
electrons traverse the undulator x ranges from 0 to 1 . We can now define the electron 
phase velocity as 



v(t) = = L[(k + /c 0 )P, - k] 

O I 



(3) 



The phase velocity measures the rate of change of the electron’s phase relationship 
(Q to the optical field. 
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In order that the rate of transfer of energy from the electrons to the radiation 
field is maximized, the wiggling electron and the electric field must be nearly resonant. 
Looking at Figure 1, we can explain this as follows: when the electric field is pointed 
in the x-direction, the electron trajectory needs to have a component in the x-direction. 
When the electric field phase increases by n and is pointed in the (-x)-direction, the 
electron trajectory must also have a component in the (-x)-direction. This will be the 
case if exactly one wavelength of light passes over the electron every time the 
electron traverses exactly one wavelength of the undulator. This is the "resonance 
condition." If one wavelength of light passes over the electron for every undulator 
period the electron traverses, the electron will maintain the same phase relationship 
with the light. Resonance thus implies that the phase velocity v(x), must be zero. 

Invoking the resonance condition, we set v(x) = 0 in (3) and solve for X with the 
result X = ?i 0 (1 - P 7 )/p 7 . From the definition of y, we find y 2 = 1 - p 2 - pf. In the 
helical undulator (as determined later), pf = K^/y 2 . Solving for p z we get 

P, - i - : 4-r- + • • • H) 

2r 

for y » 1. Substituting (4) into the expression for X above we find 

X = (5) 

2f 

which is the wavelength of light emitted by the electrons in the FEL. This expression 
shows that we can change the wavelength by changing the electron beam energy 
y me 2 or the undulator characteristics ^ or K. 

In the proceeding discussion we defined and then discussed resonance for a 
single electron. In a real electron beam, there are about 10 10 electrons in the 
micropulse; all of them could not possibly be exactly resonant. Although we define 

o o 

resonance as v = £ = 0, where ( ) represents differentiation with respect to x. we will 

find that in the low gain FEL (section E) it is desirable to operate at just above 
resonance v(0) = v 0 = 2.6. We can understand this by considering that about 10 6 
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electrons are uniformly spread over the length of an optical wavelength. If they are 
resonant, they maintain the same phase relationship with the light throughout the 
interaction even though the light is passing over them. Thus at resonance the 
radiation from the i th electron at £,■ will exactly cancel that from an electron at a phase 
(C/ - rc). This cancellation will occur for each electron. To avoid this cancellation, for 
net radiation, the electrons must radiate together, or "bunch" on the scale of the 

O oo 

optical wavelength. This requires that \ > 0 and £ > 0. We will discuss this 
microscopic bunching in part three of the next section. 

C. ELECTRON DYNAMICS 

Electrons in the undulator interact with both the undulator field and the radiation 
field which has been set up in the cavity. In subsection 1, we find the electron 
trajectories in the undulator using the Lorentz force equation of motion. We then use 
these trajectories in the second subsection to deduce measures of electron beam 
quality. In the third subsection we assume that electron radiation has set up a 
classical radiation field in the optical cavity and then allow it to interact with the 
electrons. This will result in the pendulum equation, which describes bunching. 

1. Electron Trajectories in the Undulator 

A linear undulator polarized in the y-direction causes electrons to wiggle in the 
xz-plane. In addition, since electron beams have a small radius and angular spread 
(which depends on the quality of the beam from the accelerator), electrons off axis in 
the y*direction are in stronger fields and are therefore focussed back toward the axis. 
This results in so-called betatron oscillations. 

The magnetic field of a linear undulator with no external focussing in the x- 
direction is given by 

S = [B x , B y S z ] = B 0 [0, sin(/c 0 z)cosh(/c 0 y), cos(/c 0 z)sinh(k 0 y)] (6) 
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where B 0 is the amplitude of the undulator field (Ref. 3). Consider an idealized 
electron which is injected into the undulator perfectly (i.e., x(0) = p x (0) = 0) in the x- 
direction. The electron motion is governed by the relativistic Lorentz force equation of 
motion in the undulator field (Ref. 5): 

(7) 



where e = \e I is the magnitude of the electron charge, and the electron’s momentum 
is p = ymc$, m is the electron mass, and c is the speed of light in vacuum. Note that 
cgs units are used. Substituting (6) into (7), noting that dy/dt = 0 since the magnetic 
field can do no work on the electron, we get the set of equations 



P* = 



eB o 
y me 



[p z sin(/c 0 -z(0)cosh(/c 0 y(0) - p y cos(/c 0 z(f))sinh(/c 0 y(0)] 




eB Q 
Y me 



p x cos(/c 0 z(0)sinh(/c 0 y(0) 



Pz =- 



eB 0 
y me 



p x sin(/c 0 z(f))cosh(/c 0 y(f)) 



(8a) 

(8b) 

(8c) 



where ( ‘ ) = d/dt. The right side of (8a) is an exact time derivative and can be 
integrated immediately with the result: 

P x = -~^-cos(k 0 z)cosh(k 0 y) (9a) 

where K = eBkQ,'2nmc 2 and B = Bq/^2 is the root-mean-square (rms) magnetic field 
strength. The dimensionless parameter K is the undulator parameter and is an 
important quantity in FEL physics, as shown in (5). The constant of integration has 
been set to zero indicating our assumption of perfect injection in this direction. The 
time dependence of the dynamical variables has been suppressed for brevity. 
Substituting (9a) into (8b) and (8c) we get 



Py 



cK 2 k 0 

~r 



cos 2 (k 0 z)sinh(2/( 0 y) 



(9b) 
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cK 2 k 

P* = — j-^sin(2/c 0 z)cosh 2 (/c 0 y) 

T 



(9c) 



Equations (9) are coupled, nonlinear and too difficult to solve analytically. In 
order to investigate the xz-plane motion we will now add the restriction of perfect 
injection in the y-direction also. That leaves us with 



Px= _^cos(, 0 z) 



p„(0) = y(0) = o 



Pz = 



cK 2 k 



Y 2 



0 sin(2/c 0 z) . 



(10a) 

(10b) 

(10c) 



Now in the FEL, K = 1 and y» 1, so that K/y 1 (Ref. 3). This allows us to expand 
equations (10) for the trajectories, keeping only terms to order (K7y) 2 . To zeroth order 
in (K7y), z(f) = cj5 z f where we take to be the average electron velocity in an 
undulator period. Using (4), we find c$ z = c(1 - (1+K' 2 )/2y 2 + ...) = c. Substituting 
z(f) = cjS 2 f into (10c) to find the next higher order contribution we get 



z(t) = cp 2 f - 



K 2 



4/foY 2 



sin(2/c 0 cP z f) + 



( 11 ) 



which is already to order (K7y) 2 . To find x(t), we will substitute the first term of (11) 
into (10a) since using the second would lead to higher order terms in x(f). We get 



Jn is 

x(t) = — — ^sin(/c 0 cP 7 f) + 
KqY 



( 12 ) 



where in both (11) and (12) we have set = 1 in the denominator. 

We can rewrite (11) and (12) in terms of the dimensionless time t (see (2)) as 



Kkq 

x(x) = — =-^-s\r\(2nNx) + • 
\2rcy 

z(t) = - — -js\n(AnNx) 

Sny 2 



(13a) 

(13b) 
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where N = L/X 0 is the number of undulator periods. As expected, the x-component of 
the trajectory is the wiggling motion caused by the undulator with one oscillation per 
undulator period. The z-component contains the average velocity term and a term 
with a fast oscillation (twice per undulator period). The oscillation is superimposed on 
the average motion because as an electron with constant speed makes an excursion 
to the side, the z-component of its velocity must slow, leading to an oscillation in z- 
trajectory. An electron makes this excursion twice per linear undulator period, which 
causes the z-trajectory to oscillate twice per undulator period. 

Now that we have found the trajectories in the x- and z-directions (13), we will 
find the motion in the y-direction. The y-trajectories are called betatron oscillations 
and are caused by the focussing action of the undulator magnetic field. We need to 
use a more realistic beam in this direction, but maintain a small enough beam such 
that k 0 y<& 1 and sinh(2k 0 y) = 2 k 0 y. If we then re-write (9b) in terms of t, the result is 

00 2 K 2 L 2 kn 

y(x) = ^ -cos 2 (2nNx)y(z) . (14) 

T 



Since a typical betatron oscillation takes several undulator periods to complete, we can 
average (14) over several periods (Ref. 3). We are then left with 

' 2 / 2 



K*L*k$ /v aVp(y) 

y(x)= — — y{x) =—§r 



(15) 



which is the equation of motion for a simple harmonic oscillator. The solution, and 
thus the y-trajectory, is sinusoidal. The harmonic potential is Vp(y) = c ojfy 2 ^ where 
cop = KLk 0 'y = 2nNKiy is the dimensionless betatron frequency. In a typical undulator 
N = 100. K - 1. y = 100. and we get o>p = 2k, or one betatron oscillation along the 
undulator. 



We have only treated betatron motion in one direction, but there can be 
betatron focussing in the other direction as well. This could be provided by external 
focussing magnets, or by machining parabolic pole-faces on the undulator magnets 



13 



(Ref. 3). Neither the amount of focussing nor the electron beam quality have to be 
the same in both directions. 



2. Electron Beam Quality and Emittance 

By inspection of (15) and knowledge of Hamiltonian mechanics, we know the 
harmonic betatron oscillations have a constant of motion 

o o 

Hp = [copy 2 (x) + y 2 (x)]/2 = [copy 2 (0) + y 2 (0)]/2. It can also be shown that when an 

o 

electron is injected off axis (y(0) = y 0 ) or with an initial angle (y(0) = LQ y where 0 y is 
the angle from the z-axis), its z-velocity is decreased by A|3 Z = -Hp/L 2 relative to a 
perfectly injected electron (Ref. 3). Taking differentials of (3) we find Av = Z./cA(3 z since 
k » k 0 . Therefore, an imperfectly injected electron experiences a decrease in phase 
velocity given by 

Av = + ^02) (1 6) 
L 1 + K* 



where we have used (5). 

In a real electron beam from an RF LINAC there can be =10 10 electrons per 
micropulse, each of which has a phase velocity change away from the average given 
by (16). Numerically we don’t simulate all 10 10 electrons, but instead distribute many 
electrons in both injection position and angle. They are distributed as a Gaussian, 
which has the distribution function 

f «,). exp(-qW) (17) 

V2rt a 



where o is the standard deviation and q is the variable being distributed (Ref. 3). 
Using (16) as a guide we define the standard deviation of the angular distnbution as 



_ 2JTA//0 2 

6 1 + K 2 



(18) 



where 0 is the root mean square (rms) initial angle from the z-axis. The standard 
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deviation of the position distribution, also called the dimensionless electron beam 
radius, is 







(19) 



where T 0 is the rms initial distance from the z-axis (Ref. 3). The significance of the 
factor VXL In, the characteristic optical mode size, will be discussed in section G of this 
chapter. From (16) we see that 



^ve 



2hNK 2 I<qTo 

1 + K 2 



( 20 ) 



is the standard deviation of an exponential spread in phase velocity caused by the 
initial Gaussian spread in position characterized by (19). The Gaussian spread in 
injection angles characterized by (18) also leads to an exponential distribution in 
phase velocity with the same standard deviation (Ref. 3). The exponential distribution 
arises because the phase velocity of an electron decreases no matter what the 
direction of the (Gaussian distributed) initial position or angle. The distribution function 
for the exponential is given by 

,so (21) 

G 

where q is the variable being distributed, i.e., v,, the phase velocity of the i th electron 
(Ref. 3). 

The injection angles and positions of the electrons in the beam are governed 
by a quantity called the "emittance." which is a measure of the quality of the electron 
beam. The rms emittance is defined as = Jtr o 0, and is a fixed quantity throughout 
the entire length of the undulator (Ref. 2). That is, one cannot change the radius of 
the electron beam without affecting its angular spread and vice versa. Solving (18) 
and (19) for 6 and r 0 and substituting these into the definition of emittance yields 

E rms (22) 
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where we have defined the dimensionless emittance e as 

e = o e Vc^ . (23) 

This dimensionless emittance will be used extensively in Chapters IV and V. 



3. The Pendulum Equation 

We now consider the interaction of the electrons with both the undulator and 
radiation fields. We have already found the macroscopic electron trajectories in the 
undulator. Our goal in this subsection is to find an expression which describes the 
evolution of bunching. We will use the field of a helical undulator in this section as it 
makes the problem mathematically easier. At the end we will state a prescription for 
transforming the pendulum equation to that for a linear undulator. 

The magnetic field on the axis of a helical undulator is given by 

4 = B 0 [cos(k 0 z), sin(/c 0 z), 0] (24) 

and the radiation field present in the cavity can be represented by 



E: r = E(x,y,z,t)[ cosy, -siny, 0] 


(25a) 


£? r = E(x,y,z,f)[sin\j/, cosy, 0] 


(25b) 



where E(x, y, z, t) is the wave amplitude and y = kz - tof + p (Ref. 3). 

To describe this problem of the interaction of electrons with the undulator and 
radiation fields, there are five coupled equations: 



if = - e (£ + ISxB) 

at 



(26) 



Si = ~e-j? 

dt me 



Y = 



1 






(27) 

(28) 



but only four unknowns: p x , p yt P 7 , and y, so we can neglect the z-component of (26). 
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The fields £ and E$ in (26) and (27) are the sum of the radiation and undulator fields. 
Substituting (24) and equations (25) into the x- and y-components of (26), we get 



d_ 

dt 



(7$l) = - 



e 

me 



£(1 - p z )[cos\|/, — sinvj/, 0] + p z e 0 [-sin(/c 0 z), cos(/c 0 z), 0] 



(29) 



where we have used p x = y and E = E(x,y,z,t). The right side of (29) contains 
two terms, the first of which is proportional to (1 - p z ) and is surely small for relativistic 
electrons. We will neglect this term and keep only the dominant undulator term (i.e., 
that which is proportional to B 0 ). We can then integrate (29) and set the constant of 
integration to zero indicating perfect injection. The result is 



= — ~[cos(/c 0 z), sin(/c 0 z), 0] (30) 

where K = eB 0 Xo/2nmc 2 is the undulator parameter since B = B 0 in the helical 
undulator. 

As mentioned earlier, £ is transverse and couples with only the transverse 
components of $, that is, (30). Therefore substituting (30) into (27) we are left with 



Y - 



eEK 

ymc 



cos(C, + <))) 



(31) 



where £ is given by (1) and <J> is the phase of the optical field as discussed above. 

As we noted earlier, (28) can be arranged as y -2 = 1 - p 2 - pf. If we 
substitute (30) for p 2 we get 




1 - Pi 

1 + K 2 



(32) 



which shows that a change in y does not result in a change in p a . only in P 7 . Taking 
the time derivative of both sides of (32) results in 



i = 

y 1 + K 2 



( 33 ) 



Construction of y'yfrom (31) and equating that with (33) results in 
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( 34 ) 



v 1 + K 2 
^ " T 2 



(/f + /c 0 )c 

Pz 



^cos(C + (1)) 
■fmc 



where we have used $ z = £/(/c + /c 0 )c. We will now make approximations using 
k » /c 0 and (3 Z = 1 , and substitute for k using (5) to get 



•• 2k 0 eEK 

5 = — \ cos(t; + <t>) . (35) 

■ft v 

We can make (35) dimensionless by setting d 2 /dt 2 = (clL) 2 d 2 ldx 2 , which leaves us 
with 

oo 

C, = la lcos(£ + <J>) (36) 



which is the pendulum equation. We have defined the dimensionless optical field 
amplitude la I as 



, , 4nNeKLE 
I 3 I — — — 

■fmc 2 



(37) 



The pendulum equation appropriate for the linear undulator can be found by the 
substitution K KV 0 (£) - ^i(£)] m (37) (Ref. 3). Here J 0 (£) and J^) are Bessel 
functions of the first kind of order 0 and 1 and £ = K 2 / 2(1 + K 2 ). 

oo 

Equation (36) shows that electrons with phase -n/2 < (£ + <f>) < nJ2 have ^ > 0 

OO 

while those with tc/2 < (£ + <f>) < 3nJ2 have £ < 0. This difference in phase 
acceleration shows how electrons which are initially distributed randomly from 
(£ + <!>) = —tc/ 2 to (£ + q) = 3nJ2 are driven together, i.e., bunched, near (£ + $) = jc/ 2 by 
the FEL interaction. In the next section, we will see what affect these bunched 
electrons have on the optical field. 



D. THE FEL OPTICAL WAVE EQUATION 

Knowing through the pendulum equation that the electrons will be bunched by 
the FEL interaction, we now find what affect those bunched electrons have on the 
classical radiation field. We use Maxwell's equations driven by the electron current to 
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determine the form of the FEL optical wave equation (Ref. 3). 

Electrons in the undulator emit spontaneous radiation because they are made 
to accelerate. After approximately one pass, there is enough radiation in the optical 
cavity to take on a classical form (Ref. 3). Over the next several passes this 
spontaneous emission spectrum is narrowed. In many cases, the final FEL emission 
spectrum, or linewidth, is determined by the length of its optical pulses through the 
Fourier transform (Ref. 3). Short optical pulses result in broader spectra than long 
optical pulses. At this point in the operation of the FEL, the radiation wavelength is 
given by (5). Since we know that the emission spectrum of the FEL is narrow (it is a 
coherent light source), it must be that the amplitude and phase of the radiation field in 
the cavity are slowly varying in both time and space. We will invoke the slowly-varying 
amplitude and phase approximation in the following derivation of the FEL optical wave 
equation (Ref. 3). 

The fields in the cavity of an FEL with a helical undulator are given by 
equations (25). Since the amplitude and phase are slowly varying in time, we can 
write the vector potential for these fields as 

X = e '» £ (38) 

(0 



where i = [-/, 1,0] is the polarization vector and \y = kz - (at + <f». The potential for 
the physical fields is given by the real part of (38). The inhomogeneous wave 
equation for is given by 



V 



2 



_Lii 

c 2 at 2 

/ 




(39) 



where J is the driving current density (Ref. 5). The slowly-varying amplitude and 
phase approximation can be written as inequalities: E « kE, p' « fc<t>. E « oE and 
<> « too with ( )’ = 9( )/Dz (Ref. 3). In the wave equation, this approximation has the 
effect of favoring first derivatives over all other terms. Substituting (38) into (39) and 
carrying out the z- and t- derivatives we are left with 
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( 40 ) 



V, 2 + 2 ik 


j_ i a] 


4- 


dz c dt 

^ J - 



E(x,y,z,t)e^ x -y ’ z ■'> = -^-7-t*e-^ kz ~ 

c 2 



where V 2 = d 2 ldx 2 + d 2 /dy 2 and t* is the complex conjugate of £. 

The driving current 7 is given by 

7(x,y,z) = -ec £ 5(x - x,) 5 (y - Vi ) 8(z - z,) ft = -ec £ 5< 3 >(T - T,) ft 



(41) 



where 5(x) is the Dirac delta function, T is the position vector of the field point 
(x, y, z), ~fj is the position vector of the i th electron, and the index / runs over all 
electrons in the micropulse. Substituting (30) and (41) into (40) and dividing by 2 ik 
results in 



(.Ay 2 + 


_a_ i_a_l 


4 + 


dz C dt 



Ee 1 * = -2n eK £ 






Y/ 



S< 3) (T - 7?) 



(42) 



where we have used (1) and E = E(x,y,z,t). 



We made the left side of (42) slowly varying by invoking the slowly varying 
amplitude and phase approximation. The right side is not slowly varying: it contains 
many Fourier components due to the presence of the delta functions which are each 
infinitely broadband. This means that while the left side describes an optical field that 
varies slowly, the right side can describe an optical field which varies arbitrarily fast. 
To smooth the right side we must replace the sum over all electrons by an average 
over some sample electrons. Since the fields in the cavity are slowly varying over 
several wavelengths by assumption, we only need to sample the electrons in one 
section of the micropulse one wavelength of light long and use the local charge 
density to weight the average. The result is 



-1-Vf + 


d_ + 1 a_V 


4 n 1 


dz C dt 



Ee'* = -2 neK 



PU 



- cp 7 f) 



< 




U - cp,t) 



(43) 



where p(z - cP 7 f) is the charge density of the traveling electron pulse, and < ••• > is 
an average over sampled electrons. 
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We can write (43) in a much simpler form if we change variables to 2 = z - ct 
and use the dimensionless time x. The coordinate 2 then follows a point on the 
optical pulse traveling at speed c from which the electron pulse slips back at speed 
c( 1 - p z ). In this coordinate system (let z 2) we find that L(d/dz + Mcd/dt) = d/dx. 
We can make (43) dimensionless by making the replacement x,y -» Vi/XT x,y where 
'IXL/n is the characteristic optical mode size. The result is 



“Vf + 

4 + 8x 



-'C 



Ee = -2neKL p < - — > 



(44) 



where we have assumed a long electron pulse so that p can be considered constant. 
Also, if the FEL efficiency is small so that changes in y are small, (44) takes the form 




a = -j<e l! >> 



(45) 



where a = la le'^, la I is given by (37), and the dimensionless current density j is 



. = 8A/(e7tKl) 2 p 
^ y’mc 2 



(46) 



To get the wave equation for an FEL with a linear undulator, we make the same 
substitution as before, K — > K(J 0 ( - J-\(£,)), in both (37) and (46) (Ref. 3). 

Now we will consider the physics contained in (45). First, we see that the 
growth of the field amplitude 5lal/cfr depends on both the current j and the wavefront 
curvature through the diffraction operator Y^. This unnecessarily complicates our 
discussion, so we will assume for the moment that the light is constructed of plane 
waves so that E(x,y,z,t) -> E(z.f) and Yfa -» 0. Then (45) takes the simple form 

o r 

a = - j<e~'"> (47) 

o 

We see that if j - 0, then a = 0 and the field does not grow. Also, if there is no 
bunching, the electrons are randomly distributed throughout each wavelength of light, 
<exp(-/^)> = 0 and again the field does not grow. Equating the real and imaginary 
parts of (47) leads to 
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(48a) 



-fr- = -y<cos(C + <w> 

la 1-^- = y'<sin(C, + <)>)> (48b) 

ox 

which shows that the field amplitude la I grows when the electrons are bunched near 
(C, + <(>) = n whereas the phase <|) is driven when bunching is near (C, + <|>) = k/2. 

The wave equation (45) and the pendulum equation (36), when coupled, allow 
the field and the electron phase space positions to evolve self-consistently. The 
pendulum equation describes the bunching of electrons in the presence of the 
undulator and radiation fields, while the wave equation shows how much the radiating 
bunch changes the optical field. The two equations form a feedback loop which under 
the right starting conditions leads to gain, our next topic. 

E. GAIN 

Now that we know how the electrons affect the optical field and vice versa we 
can translate this interaction to gain. FELs can be separated into two regions: low 
current/low gain defined by ; <sc 1 and high current/high gain defined by » 1 (Ref. 3). 
As it turns out, the low gain behavior we will analyze has a broader range of 
applicability, up to at least j - 1. An FEL with j = 100 will exhibit high gain 
characteristics. The difference between the low and high gain regions is in the growth 
rate of the optical field a (see (47)). We will find that the gain mechanisms in these 
two regions are different and therefore warrant a distinction between the two. 

Gain for either high or low current is defined as 

q - — 1 ~ a ° _ Energy Lost by Electrons 
Optical Field Energy 

where the last equality anses from conservation of energy. Here a 0 is the optical field 
strength seen by the electrons upon entering the undulator. Using (3), (4), and (5) it is 
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not difficult to show that the change in electron phase velocity due to a change in 
energy Ay/y or wavelength emitted AX/X is 



Av = AkN-^X- , 

y 



(50a) 



or 



Av = 2nN-— , 

A 



(50b) 



so that the energy lost by one electron can be expressed as 




(51) 



where E = y me 2 is the electron energy. The number of electrons in a volume of the 
optical mode dV is pFdV, where the filling factor F is the ratio of the cross sectional 
areas of the electron beam and optical mode. The optical energy contained in that 
same volume dV is E 2 dV/4n (Ref. 5). After some algebra we can write (49), letting 
j -> jF, as 



which is valid for both low and high current. An average over the phase velocities v, 
in the sampled portion of the electron beam is indicated by < ••• >. 

In the two subsections which follow we will discuss gain in the low current and 
high current FELs, respectively. We will find expressions for the single pass gain in 
both cases as a function of x. In FEL design, low gain usually implies an oscillator 
configuration while high gam (e.g. j = 10 4 ) implies an amplifier in which electrons 
execute only one pass before the light is outcoupled. In this case, the dimensionless 
power P = la' 2 evolves as P = (1 + G)P 0 and a large G ensures P » P 0 . After n 
passes in the oscillator P n = (1 + G) n P 0 so although G is small, over many passes 
the light is amplified until P n » P 0 . 




(52) 
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In the following analysis the wave equation is in the form (47). In section G we 
will discuss the effects of diffraction on the FEL interaction. 



1. The Low Gain FEL 



In exact analogy with a simple pendulum, the tool used to analyze the slow 
electron dynamics (bunching) in the FEL is their (£,v) phase space, one 2k section of 
which is one wavelength of light long. 



In the low gain FEL, j « 1 implies through (47) that changes in the field are 

o 

small. Multiplying both sides of (36) by £ and rearranging terms leads to 



_d_ 

dx 



— - la lsin(C, + q) 



= 0 



(53) 



and we recall that v = This has a simple solution 



v(C)=± 



2 H q + 2la lsin(£ + 0) 



1/2 



(54) 



where H 0 = Vq/ 2 - la lsin(C, 0 + 0) is a constant of the motion for each electron starting 
at (Co- v o) ar, d determines which phase space path v(Q it will traverse. These paths 
are plotted in Figure 2. For electrons starting on resonance v 0 = 0, (54) predicts fixed 
points at (£ 0 + <(>) = k/2 and 3n/2. Analysis shows that the fixed point at n/2 is stable 
and surrounded nearby in phase space by harmonic orbits (Ref. 6). .The fixed point at 
3kJ2 is unstable and defines the path characterized by H 0 = la I, called the separatrix, 
which is 



v s = ±\'2la : 



1 + sin(£ + <>) 



1/2 



(55) 



The separatnx separates phase space into regions of closed and open orbits and has 
a peak to peak height of 4>.7aT (Ref. 6). It is plotted in bold in Figure 2. 



In Figure 3 we have distributed 20 resonant (i.e., v 0 = 0) electrons uniformly 
through phase space and show their evolution as x : 0 1 . At t = 0, the electron 

position is gray and steadily darkens as t -♦ 1. The phase space simulation solves the 
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wave and pendulum equations numerically with j = 1 and initial field amplitude a 0 = 3. 



FEL PHASE SPACE PATHS 
a=l <J)=0 




Figure 2. The phase space paths v(Q traversed by the 
electrons. The separatrix is drawn bold. 



Since the electrons are resonant, they start and remain inside the separatrix on closed 
orbits. One can clearly see in Figure 3 that the electrons become bunched at £ = nf 2, 
as predicted at the end of section C. Even though the electrons are bunched, 
however, gain is nearly zero because the electrons are bunched at the wrong phase 
for field amplitude growth We can see in the figure that as many electrons move up 
in phase space and gam energy from the light wave (see (50a)) as move down in 
phase space and lose energy to the light wave. Something must be done to upset this 
balance in order to get a net gain. 

It is interesting to note that Figure 3 shows a very small amount of gam. This 
is due to the slight phase shift given by (48b) which causes the bunch to move just 
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slightly toward (£ + <(>) = n. 



*** FEL Phase Space Evolution *** 



j=l a o“3 v o =0 




Lmuumiiiiiimiiiimimumui 



-71/2 £ 371/2 0 T 1 



0.13 



0 



0.1 



0 



Figure 3. Phase space simulation of low gain FEL with 
electrons starting on resonance. Gain is negligible. 



To get an analytic solution for the gain in a low current FEL, expand (36) in 
weak fields la I = a 0 n with the result 



v(t) = v 0 + — 



>0 L 



sin^o + v 0 x) - sin(; 0 ) 



(56) 



a lL 1 

v 0 3 1" 4 



cos(2; 0 + 2v 0 t) - cos(2^o) + cos(v 0 t) - 1 - v 0 tsin(; 0 )cos(; 0 + v 0 t) 



where a 0 = la (0) I and 0(0) = 0 (Ref. 3). When we perform the average in (52) only 
the term proportional to a l in (56) survives and the result is 

2 - 2cos(v 0 t) - v 0 Tsm(v 0 x) 



G(x) = y 



vo 3 



(57) 
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We see that the gain at the end of the undulator (x = 1) depends only on the initial 
phase velocity and current density when the system is in weak optical fields. 

In Figure 4 is the result of a simulation which solves the pendulum and wave 
equations numerically for many different initial phase velocities and plots the final gain 
G(x = 1) and phase shift <))(x = 1) versus v 0 . The v 0 axis can be viewed as either a 
function of electron energy or emitted wavelength through equations (50). In the 
figure, v 0 ranges from -12 to 12, j = 1, and a 0 = 0.001 <?; n. The gain spectrum is 
anti-symmetric in agreement with (57) and has a maximum of G = 0.13 at v 0 = 2.6. 
The phase curve shows that the shift is maximum for resonant electrons, as discussed 
above, and is quite small where the gain is large. 



**** Gain and Phase Curves *** 



j=l 



a =0.001 
o 



0.13 




-12 



- 0.08 



12 



Figure 4 Gain and phase spectra for low gain FEL. 
Gain is maximum at v 0 = 2.6. 
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We have in Figure 5 another phase space simulation using the same 
parameters as in Figure 3 except that v 0 = 2.6 in order to maximize gain. In this case 
only about half of the electrons are trapped in the separatrix while the other half are 
situated on open orbits. A clear bunch of electrons can be seen at C, - n which affords 
maximum gain with a small phase shift <J>. Compared with Figure 3, more electrons 
lost energy to the light wave than absorbed energy from the light wave. It is also true 
that if the electrons started at v 0 = -2.6 the result would be maximum absorption as 
more electrons would be absorbing energy than losing energy. This behavior could 
have been predicted from studying Figure 2: since all electrons travel nearly the same 
distance along their orbit (from (36)), starting them higher in phase space (i.e., v 0 > 0) 
results in electrons with C, 0 < tc/ 2 traveling more or less parallel to the C, axis, while 
those with £ 0 > rc/2 travel down, losing energy. 




-71/2 C, 



3 rt/2 0 T 1 



Figure 5 Phase space simulation of low gain FEL with 
electrons starting at v 0 = 2.6 for nearly maximum gain of 
G = 0.12 with a 0 = 3. 
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2. The High Gain FEL 

In the high gain FEL, / » 1 and changes in the optical field are not small. 
Energy conservation cannot be used to determine the gain in this case because there 
is a substantial change in both the optical field amplitude and phase. It is still possible 
to obtain an analytic expression for the gain in weak optical fields; we merely state the 
result: 

G(t) = |-exp[(l/) 1/3 V3'c] (58) 

where v 0 is set to zero (Ref. 3). Unlike the low gain FEL (57), here the gain is 
exponential in x. We can ascertain the dependence on Vq by using the same gain and 
phase simulation used in Figure 4. Shown in Figure 6 is that simulation with j = 100 
and a 0 = 0.001 <^c n. 



**** Gain and Phase Curves ***+ 
j=100 a o =0.001 



71 

- 0.0 
j-71 




-12 v 0 12 

Figure 6 Gain and phase spectra for high gain FEL. 
Gam is maximum at v 0 = 1.4, but is still substantial at 
v 0 = 0 where the phase shift <j> = rc/2. 
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The maximum gain is G = 71 at v 0 = 1.4, closer to resonance and much higher than in 
the low gain case (Figure 4). Note that the phase shift at maximum gain is quite large. 
The gain spectrum is nearly symmetric, and substantial gain exists at resonance 
where the phase shift 4> = n! 2. 

To understand the mechanism for high gain, a phase space simulation using 
y' = 100, a 0 = 1, and v 0 = 0 is shown in Figure 7. The electrons which start at 
resonance drive the optical phase <)>, represented in phase space by the shifting 
separatrix, and attempt to bunch at £ = n/2 as discussed earlier when considering 
Figure 3. 



*** FEL Phase Space Evolution *** 

j = 100 a =1 v =0 

o o 




-71/2 C, 371/2 0 T 1 

Figure 7 Phase space simulation of high gain FEL with 
electrons starting at resonance. Final gain G = 62 = e 4 
and the phase shift o = tl'2. The separatrix is shown 
growing with the field strength la I and shifting left as the 
phase o grows. 
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It is the phase shift <{> = kJ2 which allows gain as the relative electron bunch phase 
(C + <>) = 7i drives the optical field amplitude. This leads to exponential gain (linear on 
the logarithmic plot of Figure 7) as predicted by (58). 

To summarize, when y <k 1 we must start the electrons at v 0 > 0 in order that 
they bunch near (C, + <[>) = n and gain results. For an FEL with /' » 1, electrons at 
v 0 = 0 drive the optical phase with the high current until the bunch is situated near 
(£ + <!>) = 7i. In the case of low gain, the optical phase shift is unimportant and small; in 
the case of high gain it is a crucial aspect of the interaction. 

3. Gain Degradation Due to Electron Beam Quality 

In section C.2 of this chapter, we discussed transverse electron beam quality 
and deduced two measures, o 0 and o ve , from the emittance. The longitudinal beam 
quality has a measure 

o G = 4k N ( 59 ) 

Y 

which is the standard deviation of a Gaussian spread in energy (see (17)), and Ay me 2 
is the rms energy difference from the resonant energy (Ref. 3). 

Consider the low gain spectrum in Figure 3 again. The full width of the main 
positive gain portion at half the maximum is approximately n and the distance from 
maximum gain to maximum absorption is approximately 2n. Therefore, if the spread in 
phase velocities of the electron beam Av = n, we should expect serious gain 
degradation as many electrons will be injected into the undulator with phase velocity 
corresponding to less than maximum or even to negative gain. Inspection of (50a) 
and (16) shows that Av has contributions from o e , o ve , and o G . 

In Figure 8 is a phase space simulation using the same parameters as those 
used in Figure 5 except that a small energy spread o G = 1 has been included. Only 
the final phase space positions of the electrons have been plotted. The electrons are 
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spread through phase space and bunching has been degraded. Gain has been 
reduced to G =0.10. 




-71/ 2 C 37T/2 0 T 1 



Figure 8. Phase space simulation of low gain FEL with 
energy spread o G = 1. Gain is reduced to G=0.10. 
Compare with Figure 5. 



In the high gain spectrum of Figure 6 we note that the gain spectrum bandwidth 
for y = 100 is approximately 2k. In general, this width is proportional to y' 1 6 in the high 
gam region (Ret. 3). Since there are no regions of absorption and the spectrum is 
about twice as wide as that of the low gain example, we would expect the high gain 
FEL to be more resistant to beam quality effects. However, when increasing the 
undulator length, beam quality measures o G . (59), o e . (18), and o ve . (20) increase 
proportional to N. This means that the quality of the electron beam is magnified by a 
long undulator. Since y » IN 3 . where / is the peak current, high gain is most easily 
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and often obtained by using a long undulator. In addition, higher peak current / from 
the accelerator leads to poorer beam quality, so that either way the designer tries to 
maximize gain the problem with beam quality gets worse. When trying to get all one’s 
light amplification in a single pass, little gain degradation can be tolerated. Indeed this 
may become one of the major drawbacks of using an FEL amplifier as a high power 
weapon. 

In Figure 9 is a phase space simulation using the same parameters as those 
used in Figure 7 except in this case o G = 1 . Only the final separatrix and electron 
phase space positions have been plotted. The electrons are spread throughout phase 
space and bunching has been degraded. Gain has been reduced from G = 62 to 
G = 50. 



*** FEL Phase Space Evolution *** 



j=100 a =1 V =0 0^=1 

O O G 




-n/2 C 371/2 0 T 1 



Figure 9 Phase space simulation of high gam FEL with 
energy spread o G = 1. Gam is reduced to G = 50. 
Compare with Figure 7, where G = 62. 
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In both Figures 8 and 9 we have only included the energy spread a G 
contribution to beam quality. Clearly the degradation would be even worse if we also 
included the effects of emittance. 

F. SHORT PULSES 

It was mentioned during our derivation of the wave equation that the electrons 
slip back relative to the light pulse at speed c( 1 - p z ). From the definition of 
resonance, we know that the distance between the leading edges of the electron and 
light pulses at x = 1 is NX. This distance is called the "slippage length." If the 
electron pulse length is on the order of the slippage length, the pulse is considered 
short (Ref. 3). In this case we cannot make the long pulse assumption in (44) 
because each section of the light pulse sees the charge density changing as the 
electron pulse slips past. 

In an optical resonator of length S, it takes time 2 Sic for the optical pulse to 
make one round trip. Let the desynchronism d be the distance between the light and 
electron pulses at x = 0 on each pass (Ref. 3). We normalize d and the electron 
pulse length o z - l e i(NX) to the slippage length. At exact synchronism, unbunched 
electrons slip back past the light pulse as x : 0 -» 1, bunching occurs, but only the rear 
of the light pulse gets amplified. On successive passes the delay in bunching 
becomes longer and the amplification moves further back on the light pulse. Since 
only the rear of the light pulse is being amplified while the losses (as determined by 
the resonator quality factor Q) are the same at all points in the pulse, its centroid 
moves at a speed slower than c. This is called "lethargy," and ultimately will lead the 
pulse out of the region of gain (that is. it won't overlap the electron pulse) and the 
pulse amplitude will vanish (Ref. 3). 

To compensate for the effects of lethargy a small amount of desynchronism 
must be added. This is accomplished by reducing the path length S that the light 
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pulse must travel by adjusting the length of the optical cavity. In this case, 
d = -AS/(NX) (Ref. 3). 

Simulations of an FEL using picosecond electron pulses will be discussed in 
Chapter III. There we will explore the effects of desynchronism on the power 
evolution. 



G. FEL WITH GAUSSIAN OPTICAL MODE 



In section E on gain we discussed the interaction of the electron beam with 
light approximated by plane waves. An optical cavity bounded by spherical mirrors 
supports Hermite-Gaussian (and other) modes, which are solutions of (45) instead of 
(47) and accurately describe the laser beam inside and outside the cavity (Ref. 7). 
The optical beam can be constructed from any number of these modes, but the one 
most desired in a weapons laser for its propagation characteristics is called the 
fundamental Gaussian mode, or just Gaussian mode. 

The usual Gaussian mode as given in Ref. 7 can be written in terms of our FEL 
dimensionless parameters. The result is 



a _ 30 --'WZo /<t> G (t) 

a G - — —e e 

w(x) 



where 



<t> G (x) = -tan 1 



X - X u 






~ tJ 

Zo + (* - * J 2 



is the phase shift of the Gaussian mode relative to that of a plane wave, 

2 



w^x) = 1 + 



X - X. 



(60) 



(61) 



(62) 



is the square of the mode radius, which is normalized to the mode waist radius w$ , 
and 
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z 0 = Wq (63) 

is the dimensionless Rayleigh length (Ref. 3). The mode waist position is i w . Radial 
dimensions are normalized to VlX/jt, the characteristic mode size. We see in (62) that 
z 0 is the distance over which the mode area doubles. In a symmetric optical cavity 
with spherical mirrors a distance S apart and radii of curvature R, the dimensionless 
Rayleigh length is given by z 0 = (2 SR - S 2 )' /2 /2L where L is again the undulator 
length (Ref. 7). The Rayleigh length z 0 and the waist position x w completely 
characterize the Gaussian mode in the optical cavity. 

The amplitude of the Gaussian mode la I is a maximum at x = i w where 
w 2 (t) = 1. It decreases away from the waist in both directions and is a minimum at 
both x = 0 and x = 1 (if i w = 0.5). At any time x, the amplitude decreases in the radial 
direction as a Gaussian, hence the name. The phase of the Gaussian mode 4> e (61) 
contains two terms. The first term provides a phase shift, called the Guoy phase shift, 
relative to the plane wave, and switches sign at the waist (Ref. 7). The second term in 
(61) decreases the phase shift for positions off-axis which creates the wavefront 
curvature necessary for focussing the beam at the waist. 

In the electron phase space, the term ao/w(x) causes the height of the 
separatrix to increase as x : 0 -» i w and then decrease as x : i w -» 1 . The effect is 
proportional to w _Vj (x) (from our knowledge of the separatrix height) and complicates 
the phase space paths. The phase of the Gaussian mode (61) is continuously shifted 
as x : 0 —> 1. The shift has the sign opposite to that caused by the FEL interaction in 
(48b). It is therefore represented in phase space by the separatrix shifting to the right 
as x : 0 — ► 1. Looking at Figure 2 again, we can predict that the electrons will have to 
start at a phase velocity v 0 > 2.6 in order to compensate for the shift and still bunch at 
the correct phase for gain. (Compare this to the case of the high gain FEL where the 
phase velocity had to be lowered in order to bunch at the correct phase.) 
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In order to observe the effects of the Gaussian mode on the FEL interaction, 
assume z 0 » 1, nearly plane waves, and expand (61) and (62) to first order in (1/z 0 ). 
The phase becomes <t> G = -(t - t w )/ z 0 , and w(x) = 1. The pendulum equation (36) in 

oo 

the presence of the Gaussian mode is C, = la G lcos(^ + <|> G ) which for z 0 » 1 becomes 






- a o e 



. —ilZn 



COS 



Co + 





(64) 



when expanded in weak fields to first order in a 0 (Ref. 3). Comparing (64) with (36) 

oo 

we see three distinct differences. First, the amplitude of the phase acceleration C, 
decreases for electrons off axis because they see a weaker field strength 
la I — > a 0 exp(-r 2 /z 0 ). Second, there is a shift in the initial phase of all electrons in 
the amount x w /z 0 . This is inconsequential since the electrons are initially distributed 
uniformly throughout phase space. The last difference is in the phase velocity, which 
is decreased by Av = 1/z 0 . This leads to a wavelength change of AXJX= V(2nNz 0 ) 
from (50b). This also requires an increase in the initial phase velocity for maximum 
gain to v™ 3 * = 2.6 + 1/z 0 in the low gain FEL. This effect is caused by the optical 
phase shift (61) and is shown in Figure 10. The figure is the result of a simulation 
similar to that in Figures 4 and 6 except that the finite size of the optical mode is 
included. Here the Rayleigh length z 0 = 5 and maximum gain is found at 
v 0 = 2.9 = (2.6 + 1/5). 

The presence of the Gaussian optical mode produces measurable effects in the 
performance of the FEL. This is the theme of the bulk of this work, including the final 
chapter on Compton backscattering where we calculate the photon emission rate for 
electrons in the Gaussian field of a laser. 
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Figure 10. Gain spectrum for low gain FEL with Rayleigh 
length z 0 = 5. Maximum gain is achieved at v 0 = 2.9. 
shifted from v 0 = 2.6 by the Gaussian mode. 
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III. SIMULATIONS OF FIREFLY 1KW FEL 



A. BACKGROUND 

The FIREFLY (Far Infrared FEL) IkW FEL at Stanford University was proposed 
in 1994 as an upgrade to the existing Infrared FEL at Stanford University (Ref. 8). 
With four new 10 MeV accelerator structures and a 5 MeV electron injector, an 
electron beam with energy 45 MeV could be generated. The electron injector 
operating at 1 1 .8 MHz would produce an average current of T = 1 ma, and the 
average power in the electron beam would be P e = 7V = 45 kW. A good measure of 
the maximum efficiency of an FEL, i.e., how much energy one can extract, is the gain 
spectrum bandwidth. Considering Figure 4 and (50b), we find that Ay/y = M(2N). If 
we assume that we can extract the full Ay/y= 1/(2/V) from the beam with N = 25, then 
the average power out of the cavity would be F = F e /(2N) = 900 W. This average 
power would be almost 100 times greater than the world’s most powerful FEL, and it 
was estimated this could be done for only one million dollars with existing technology. 

The proposed FEL and two variants with lower electron beam energies of 25 
MeV and 10 MeV were simulated. All other parameters were kept the same and are 
listed in Table 1. 

B. SIMULATIONS AND RESULTS 

Because all three FELs simulated have emission wavelengths X in the far 
infrared, the charactenstic mode size vLX/n is large enough that diffraction (i.e., the 
Gaussian mode) should not dominate the interaction. The slippage length A/?., 
compared to the electron pulse length, is also large, so the effects of short optical 
pulses must be addressed. (Some basic consideration of short pulse effects was done 
in section F of Chapter II.) To study short pulse effects, simulations following 
longitudinal multimodes were used. These simulations divide the optical pulse into a 
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number of fixed longitudinal field sites, or bins, a(z) (Ref. 3). The simulations are 
called "longitudinal multimode" because dividing the optical pulse into longitudinal 
sites a — > a(z) is equivalent to a -> a(k z ) in Fourier space; the number of sites z is 
equivalent to a number of modes k z . 





1 


II 


III 


Beam Energy 


45 MeV 


25 MeV 


10 MeV 


Peak Current 


40 A 


40 A 


40 A 


Micropulse Length 


2.1 ps 


2.1 ps 


2.1 ps 


Undulator Type 


Linear 


Linear 


Linear 


K 


1.1 


1.1 


1.1 


N 


25 


25 


25 


Xq 


6.0 cm 


6.0 cm 


6.0 cm 


X 


8.4 pm 


26.6 pm 


156.7 pm 


j 


1.7 


2.9 


7.2 




3.2 


1.0 


0.2 



Table 1. Parameters for FIREFLY IkW FEL (column I) 
and two variants (columns II and III) used in simulations. 



The slower moving electrons slip back past the optical pulse to sites z - x and interact 
with a number of field sites corresponding to the slippage distance NX (recall that z is 
normalized to the slippage distance). Conversely, each field site a(z) interacts with a 
length NX along the electron micropulse as the micropulse slips back. We generalize 
the pendulum (36) and wave (47) equations to include the slippage: 
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(65) 



C, z+x — ls z icos(^ z+T + <() z ) 

3 z = —j z _ K <e >z+x • (66) 

As x : 0 — > 1 in increments of Ax in the simulation, (65) and (66) show how the field 
changes at site z are caused by the electron dynamics at the micropulse site z + s Ax, 
where s = 0, 1, 2, ...,N t is an integer and N, = 1/Ax is the total number of time steps, 
and vice versa. 

Contained in Figure 11 is the output from a simulation of the 1 kW machine, 
whose parameters are listed in column I of Table 1. The dimensionless current 
density is j = 1.7, the micropulse length is o z = 3.2, cavity losses are determined by 
0 = 100, and the desynchronism d - 0.03 is large (see Chapter II section F for 
definition of d). 



**** FEL Pulse Evolution **** 
j=1.7 CT 2 =3.2 Q=100 d=0 . 03 




-20 z 



20 -19 v 19 0 



n 2000 



Figure 11. Longitudinal multimode simulation of 
FIREFLY. 
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The parabolic electron pulse is drawn in the lower left window at x = 0 and at x = 1 to 
illustrate slippage. Above the electron pulse is plotted the optical field amplitude 
la(z,n)l at each pass n. The scale ranges from zero in black to the maximum in 
white; a contour line is drawn at the median field amplitude. Above the field evolution 
is plotted the final optical field amplitude la(z,n = 2000)1 in black. At the bottom 
center is drawn the gain spectrum (57). Above this is the optical power spectrum 
evolution P(v,n), with the final power spectrum P(v,n = 2000) drawn on top in black. 
The power spectrum P(v,n) is the Fourier transform of the square of the optical field 
la(z,n)l 2 at each pass n. The phase velocity v in the power spectrum is used in the 
sense of (50b); that is, it represents the wavelength \ so that P(v,n) shows the 
linewidth of the radiation. The total power evolution P{n) is plotted at lower right in 
the figure. Above this is plotted the electron phase velocity distribution f(v,n). The 
distribution /(v,n) is a projection of the final electron phase space positions (£,v) onto 
the v-axis at each pass n. As such, v in this plot is to be understood in the sense of 
(50a); that is, it represents the electron energy y me 2 so that /(v,n) is interpreted as an 
energy spread. 

Having described all the plots in Figure 11, we now interpret their features. 
The large desynchronism has caused the optical pulse la(z,n)l to be artificially long 
by advancing the front edge of the pulse ahead of the electrons. The electron pulse 
amplifies only the trailing edge of the optical pulse on each pass. This leads to 
reduced coupling; the peak field strength is only la I = 25. As evidenced by the power 
spectrum P(v,n), the long, smooth optical pulse leads to a narrow linewidth. The 
electron phase velocity distribution f(v,n) does not evolve any further after the first 
few hundred passes. At saturation, the electrons reach the bottom of the separatrix in 
phase space (Ref. 3). This means that the electron distribution f (\\n) should have a 
width of about 4VT51 = 18 at saturation, the peak-to-peak height of the separatrix at 
la I = 20, the final field amplitude. In the figure, the width of /(v,n) is about 18, 
showing good agreement. The power P(n) at each pass n is obtained by integrating 



42 



over the optical pulse la(z,n)l 2 . The power evolution P{n) plotted at lower right 
shows that a steady state is reached after approximately n = 700 passes where 
P(n) = 91. 

Plotted in Figure 12 are the results of several simulations like that above. The 
steady state dimensionless power is plotted versus the desynchronism. 



1500 
1000 
500 
0 

0 0.005 0.01 0.015 0.02 0.025 0.03 

d 

Figure 12. Plot of outcoupled dimensionless power 
versus desynchronism for FIREFLY. 




The maximum on the curve at d = 0.0005 corresponds to an actual output average 
power P = 900 W. the desired output, but requires the injector to operate at 23.6 MHz, 
or double the rate discussed earlier. Evidently the short pulses lead to reduced 
coupling between the electrons and optical field, lowering FEL efficiency. Doubling the 
operating frequency of the electron injector was proposed in Ref. 8 as a follow-on to 
the 1 kW upgrade for the purpose of attaining 2 kW average power. Note that the 
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maximum in Figure 12 is rather sharp. Since the slightest change in d can produce a 
large change in the output, the FEL is less stable when operated in this region (Ret. 
3). 

The FIREFLY variant listed in column II of Table 1 has a shorter dimensionless 
electron pulse length o z than that in Figure 1 1 . The output of a simulation with 
j = 2.9, a z = 1, Q = 100, and d = 0.005 is shown in Figure 13. The features of the 
plots show the occurrence of the trapped-particle instability (Ref. 3). 



**** FEL Pulse Evolution **** 
j=2 . 9 0 z =l Q=100 d=0 . 005 





104 P (V,n) 



f (v,n) 
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Figure 13. Longitudinal multimode simulation for the 
variant of FIREFLY listed in column II of Table 1. 



This instability occurs when the field amplitude la I becomes so strong that electrons in 
phase space traverse their entire path (see Figure 2) during a single pass through the 
undulator, executing a synchrotron oscillation with frequency v s = \T a~ = 2 n (Ref. 3). 
In real space, this oscillation represents motion superimposed on the electron 
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trajectory and leads to emission into sidebands v 0 ± v s (Ref. 3). As these sidebands 
grow, the FEL power will increase and the fields will get stronger, causing more 
electrons to execute synchrotron oscillations, etc. At least four sidebands can be seen 
in the power spectrum P(v,n) in Figure 13. These sidebands modulate the optical 
pulse envelope la(z,n)l as shown, and the increased field strength causes the wide 
electron distribution f(v,n). In this case, the instability has imposed a periodic 
structure on the optical field envelope la(z,n)l and the power spectrum P(v,n) which 
results in oscillations in the total power P(n) over many passes n. 

Figure 14 plots the dimensionless power versus desynchronism from several 
simulations like that of Figure 13. 




Figure 14. Plot of dimensionless outcoupled power 
versus desynchronism for vanant of FIREFLY listed in 
column II of Table 1. 
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We see from the plot that Figure 12 with d = 0.005 lies in the transition region in the 
desynchronism curve, where the fields are strong enough to cause the trapped-particle 
instability, but where operation of the FEL is relatively stable. The peak in this curve 
at d = 0.003 corresponds to an outcoupled average power of F = 140 W. 

To further illustrate the effects of slippage, the example listed in column III 
represents FIREFLY operating at a much longer wavelength X= 156.7 pm. The 
slippage length is therefore long and leads to a small value of o z =0.2. Figure 15 
contains the result of a simulation using j = 7.2, o z = 0.2, Q = 100, and the very small 
value of d = 0.0005. 



**** FEL Pulse Evolution **** 
j=7 .2 O z = 0.2 Q=100 d=0 . 0005 




133 P(v,n) 



f (v,n) 



6000 



-4 z 4 -101 V 101 -101 V 101 



Figure 15. Longitudinal multimode simulation for variant 
of FIREFLY listed in column III of Table 1. 



The small desynchronism means that, unlike in Figure 11. the light pulse a(z,n) i has 
not been artificially lengthened. The optical pulse reaches steady state with a length 
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comparable to the electron pulse length shown in the plot at lower left. The fields are 
strong and the trapped-particle instability again leads to structure in the optical pulse. 
We should note here that the spikes to the left of the main optical pulse are only about 
one wavelength of light in length. Because this would violate the slowly-varying 
amplitude assumption in the wave equation (47), these may be meaningless. The 
short, modulated optical pulse results in a broad power spectrum P(v,n) (note the 
scale in v compared with that in Figure 13). The strong fields also result in a broad 
electron distribution. Maximum power P{n) = 900 is achieved at n ~ 1500 passes, but 
steady-state P{n) = 520 isn’t reached until n ~ 4000 passes which is predicted to take 
about 0.3 ms in FIREFLY. 

Figure 16 plots power versus desynchronism for column III of Table 1. 
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d 

Figure 16. Plot of outcoupled dimensionless power 
versus desynchronism for variant of FIREFLY listed in 
column III of Table 1. 
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We note that d = 0.0005, that used in Figure 15, is the location of the peak in this 
curve and corresponds to an average power of P = 17 W. 

C. CONCLUSIONS 

The operational characteristics of far infrared wavelength FELs driven by short, 
picosecond electron micropulses will be determined by short pulse effects. The 
desynchronism coupled with the cavity Q will affect the length and structure of the 
optical pulse. This in turn determines the width and structure of the power spectrum. 

Based on our simulations of the 1 kW FEL, as depicted in Figures 11 and 12, 
the goal of 1 kW average power can be achieved, but at double the electron injector 
duty cycle first proposed. More simulations should be done to substantiate this 
possibility. Doubling the injector duty cycle was proposed in Ref. 8, but as a follow-on 
to the 1 kW proposal as a means of achieving 2 kW. We also note that the trapped- 
particle instability will probably occur. 
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IV. TRANSVERSE EFFECTS IN UV FREE ELECTRON LASERS 



A. BACKGROUND 

In an ultraviolet free electron laser (UVFEL), the characteristic optical mode 
size VTaTtc tends to be small because of the short optical wavelength X. For example, 
if X = 200 nm and L = 3 m, 'U.X/n = 0.4 mm. This is also approximately the size of 
the electron beam. For the ultraviolet FEL, the effects of the spatial extent of the 
electron beam as compared to the optical mode are of concern. 

The UVFEL under consideration in this chapter is that proposed by the Laser 
Processing Consortium (LPC) at Thomas Jefferson National Accelerator Facility 
(TJNAF) in Newport News, Virginia (Ref. 9). The LPC is a group of industry and 
university scientists and engineers interested in an efficient, reliable, coherent UV light 
source for processing materials. On the order of 100 kW of UV light is needed to 
process commercial products ranging from nylon to steel. The parameters for the 
UVFEL proposed in Ref. 9 are listed in Table 2. 



Beam Energy 


200 MeV 


Peak Current 


270 A 


Micropulse Length 


0.25 ps 


Undulator Type 


Linear 


K 


0.95 


N 


72 


Xq 


3.3 cm 


X 


204 nm 



Table 2. Parameters for the LPC proposed UVFEL. 
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In this chapter we analyze the LPC UVFEL using transverse multimode 
simulations. These simulations are described in section B. We use these simulations 
in section C to study the effects of the optical mode size on FEL gain. In section D we 
study the effects of electron beam focussing on FEL gain. 



B. TRANSVERSE MULTIMODE SIMULATIONS 



In Chapter III we extended the pendulum and wave equations in the z-direction 
in order to properly simulate longitudinal multimodes by dividing the optical pulse into 
z sites. To simulate transverse multimodes, we make the same transition in the x- 
and y-directions. The optical field is extended to sites in x and y so that a — > a(x,y). 
Generalizing the wave equation (45) and the pendulum equation (36) gives (Ref. 3): 



V I + 4- 

4 3x 



a (x ,y ,x) = -<j exp(-/ £(x ,y ,x))> (xyiT) 



C (x.y,*) = 



a(x,y,x) 



cos 



C(*.y.x) + <M*,y,x) 



(67) 

(68) 



The transverse dimensions are normalized to the characteristic mode size V\L In. The 
set of equations (67) and (68) are fully self-consistent in three dimensions (x,y,x). 

For the UVFEL represented in Table 2, we calculate an electron pulse length 
slightly greater than five slippage distances o z = 5.1. Therefore, it can be assumed 
that slippage will not be an important consideration, so that we may follow only one 
representative longitudinal site as in (67) and (68). 

In section G of Chapter II, we discussed the fundamental Gaussian mode but 
mentioned that in general any number of Hermite-Gaussian transverse modes are 
present. In order to model the evolution of the FEL from weak fields through 
saturation, these higher-order modes are included as the FEL interaction amplifies and 
distorts the optical mode. Both the optical field amplitude and phase can be altered 
(Ref. 3). An example of this is shown in Figure 17. 
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Figure 17. Optical multimode evolution from weak fields 
through saturation. 
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Figure 17 shows the evolution of the FEL from weak fields a 0 = 0.0001 through 
saturation using j - 20 and a small energy spread o G = 0.5. The Rayleigh length for 
the cavity is z 0 = 0.5. The top set of plots show the FEL oscillator after n - 12 
passes. The electron phase space plot shows no noticeable bunching and the final 
gain G = 0.94 is in steady state (i.e., G(n ) constant ) with power increasing 
exponentially. At right is shown the optical field amplitude between the mirrors 
la(x, t)I. The mode is clearly distorted from a fundamental Gaussian mode with the 
peak field (in black) located at x = 1 where single pass gain is the maximum. Moving 
down the set of figures, the electron phase space begins to exhibit bunching at 
n = 24. The bunch descends along the v-axis, indicating energy loss, and near 
saturation at n = 28 it starts to ascend taking energy back from the optical field. At 
saturation (n = 40), the electron bunch is clearly moving up in phase space at x = 1. 
The motion of the electron bunch through phase space, caused by the strengthening 
fields, is attended by a drop in final gain G(x=1) as the electrons are absorbing light at 
the end of the undulator. This means that on each pass, G(x) peaks at x < 1 along 
the undulator and the peak field moves back along the undulator. Scanning down the 
page through the plots of la(x,x)l, we see the peak field strength moving back along 
the undulator. At saturation, the peak gain per pass is nearly at the middle of the 
undulator and the final gain is reduced to G = 0.0034. The optical field la(x, x)l 
resembles the fundamental Gaussian mode at saturation where the gain is low. 

The evolution of the distorted optical mode to the Gaussian requires the use of 
the transverse multimode simulations described above. Multiple transverse modes are 
included in the simulation to construct the distorted optical fields of Figure 17 for the 
same reason it takes a large number of sine waves to synthesize a square pulse in 
Fourier analysis. 
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c. 



OPTICAL MODE SIZE EFFECT ON FEL GAIN 



It has been shown that the Rayleigh length which optimizes gain for a filament 
electron beam (i.e., o e =0) is z 0 = 1/VT2 = 0.3 (Refs. 3, 10). This is the Rayleigh 
length which maximizes the average filling factor F (averaged along the undulator 
length) by minimizing the volume occupied by the Gaussian optical mode around the 
electron beam. Reference 10 also discusses gain in a Gaussian mode for non-zero 
electron beam radius, but again using arguments based on the filling factor. The 
conclusion, not surprisingly, is that maximum gain is achieved when the Rayleigh 
length is chosen such that the electron beam fits inside the optical mode (recall 
z 0 = w $ ). But in the case of the UVFEL, it may be necessary to have an electron 
beam bigger than the optical mode because of the small mode size VXT/jr. In that 
case, one cannot quantify gain in terms of the filling factor, which may be greater than 
one. The three dimensional simulations described in the previous section are used to 
find the influence of the Rayleigh length on FEL gain. 

Figure 18 shows the result of a three-dimensional, transverse multimode 
simulation of the UVFEL parameters listed in Table 2. The corresponding 
dimensionless parameters are current density j = 30 with F is no longer included in 
v 0 = 4.5 is the optimum initial phase velocity, the Rayleigh length is z 0 = 0.3, 
a 0 = 0.001 is the initial optical field amplitude, the cavity loss is Q = 5, o e = 0.45 is the 
dimensionless electron beam radius, o e = 0.9 is the angular spread, the energy spread 
is c G = 1.3, and = 1.0 is the betatron frequency. The dimensionless emittance 
r = 0.43 corresponds to an actual normalized emittance e w = 11n mm-mrad 
(r N = Yfmw) The plot at top left shows the evolution of the optical field amplitude at 
the end of the undulator with each pass la(x,n)l. To the nght is plotted the final 
optical field amplitude at the nght-hand mirror la(x,y)l. The optical field between the 
mirrors la(x,x)l at y = 0 is plotted in the center; the tick marks at t = 0 and x = 1 show 
the beginning and end of the undulator. Overlaid on this plot are the positions of a 
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random sample of the electrons used in the simulation. In this figure, cr e /w 0 = 0.8 and 
the electrons near the edges of the Gaussian mode see weaker optical fields and thus 
a degraded interaction (see (60)). 



TRANSVERSE FEL MODEL 
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Figure 18. Transverse multimode simulation of UVFEL. 
Gain is degraded due to large electron beam 
radius/optical mode waist radius ratio o e /w 0 = 0.8. 



The plots along bottom show the electron phase velocity distribution at each pass 
f(v.n), the final phase space positions (£,v) of a sample of electrons, and the power 
P(n) and gam ln(1+G(n)) evolutions. The gain G(n ) is plotted at the end of the 
undulator G(t= 1) after each pass n. These plots show FEL gain in steady state (i.e.. 
G{n) approximately constant) with final gain G = 0.75 and power increasing 
exponentially. Electrons in phase space are spread out due to Oq and o 0 so that 
bunching can barely be seen at £ = n even though the limit of weak fields la I = n has 
been reached. 
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The results of many simulations like that in Figure 18 were used to plot gain 
normalized to the maximum one-dimensional gain, G 0 = O.135/F, where 
F = jtcf/(z 0 + 1/12z 0 ) is the average filling factor, versus Rayleigh length (Refs. 3, 10). 
The result is plotted in Figure 1 9. 







Figure 19. Plot of normalized gain versus Rayleigh 
length for UVFEL. Normalized gain is minimum where 
optical mode volume is minimum z 0 = 0.3. 



The first feature to note in the figure is that for the UVFEL, z 0 = 0.3 corresponds to the 
minimum normalized gain because the optical mode volume is minimized. The optical 
mode volume must be large enough to envelope the electron beam Decreasing the 
Rayleigh length below z 0 = 0.3 causes the mode to expand quickly away from the 
mode waist. This makes the mode large at the ends of the undulator. enveloping the 
electron beam, and increases normalized gain. Increasing the Rayleigh length from 
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z 0 = 0.3 makes the optical mode larger everywhere and increases gain until the mode 
gets so large, z 0 = 0.9 in this case, that electrons see weaker fields everywhere and 
gain is reduced. For a larger electron beam, the optimum value of z 0 for normalized 
gain would be larger than z 0 = 0.9. 

D. ELECTRON BEAM SIZE EFFECT ON FEL GAIN 

Now, given that the Rayleigh length of the optical cavity has been chosen, how 
is FEL gain affected by the choice of electron beam radius? Since cop = 1 .0 for this 
UVFEL, the electron’s betatron trajectories are nearly straight lines, allowing them to 
be focussed at the undulator center using external magnets. Keeping emittance 
e = o e V^ fixed, decreasing the electron beam radius at the waist increases the 
angular spread and degrades the interaction, so there must be a limit to how small 
one can focus the electron beam. In all simulations, the Rayleigh length is held fixed 
at z 0 = 0.5. 

Figure 20 contains the output of a simulation of the UVFEL electron beam with 
e = 0.43. In this example, we have reduced the electron beam radius to c e = 0.3, 
which increases angular spread to o e = 2.1. We also were required to increase the 
current density to j = 68 in order to keep peak current constant. The optimum initial 
phase velocity is v 0 = 4.3 and all other parameters are the same as those in Figure 
18. The phase space plot (i^,v) is dominated by the (exponential) angular spread 
which indicates o e has been made too small. In the plot of la ( x,t) 1 in the center, one 
can see the electrons focussed to c e = 0.3 at the center of the undulator, but the large 
angular spread cames many electrons outside the optical mode at x = 0 and x = 1 . 
This leads to gain degradation. 

Figure 21 plots the normalized gain versus electron beam radius for three 
different emittances: (a) c = 0.28, (b) £ = 0.43, and (c) e = 0.6. (The simulation in 
Figure 20 is from curve (b).) We see in the curves that decreasing the electron beam 
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radius from a e = 1 .0 increases gain as more electrons are enveloped in the optical 
mode. The beam radius can be made too small, however, causing gain to decrease 
because of the increasing angular spread. 
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Figure 20. Transverse multimode simulation of UVFEL. 
Electron beam focussed at x = 0.5 but angular spread 
causes beam to expand outside optical mode at t = 0 
and 1. 
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The optimum electron beam radius, and the peak gain at that radius, depend strongly 
on the beam emittance. With better emittance, curve (a) in Figure 21. the optimum 
electron beam radius is small. c e = 0.3, and peak gain is increased to G/G 0 = 66%. 
Increasing emittance causes the optimum a e to grow larger and the maximum to 
become more broad until (curve (c)) the maximum gain is only weakly dependent on 
a e due to the large angular spread at any radius. This maximum gain is reduced to 
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G/G 0 ~ 10% at a e =0.6. The expected emittance in the LPC UVFEL is e = 0.43, 
curve (b). We find that the optimum electron beam radius is o e = 0.4, slightly smaller 
than that proposed. 




Figure 21. Plots of normalized gain versus electron beam 
radius for three different beam emittances: (a) e = 0.28, 

(b) e = 0.43, and (c) e = 0.6. 

E. CONCLUSIONS 

The small characteristic mode size of the short wavelength UVFEL makes the 
electron beam/optical mode overlap a crucial consideration. Making the optical mode 
larger by either decreasing or increasing the Rayleigh length from z 0 = 0.3 results in 
increased normalized gam lor a given electron beam radius. This is an important 
result for UVFELs. It shows that it is not necessary to have the electron beam 
enveloped by the optical mode throughout their interaction. With Rayleigh length fixed. 
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focussing the electron beam to a smaller radius will increase gain until the angular 
spread, which must increase as the electron beam radius is decreased, degrades the 
interaction and gain is reduced. With improvements in emittance, the optimum 
electron beam radius is smaller and the peak FEL steady-state gain is significantly 
increased. We find the optimum electron beam radius for the LPC UVFEL is o e = 0.4. 
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V. FELS WITH SHORT RAYLEIGH LENGTH 



A. MOTIVATION 

The problem of high power density on optical cavity mirrors is shared by both 
UVFELs and high average power FELs, the former due to its short wavelength and 
therefore a small diffraction rate, the latter due to extremely high (= 100 MW) 
intracavity average power. This problem could be solved by designing a long optical 
cavity with z 0 = 0.3, allowing space for diffraction to make the beam spot on the 
mirrors large. Another way to make a large spot is to use a cavity with a short 
(z 0 < 0.3) Rayleigh length. When the FEL has limited space for the optical cavity (on a 
ship, for example), the second method of mirror protection is particularly attractive. 

As an example of the effect of the Rayleigh length on mirror power density, 
consider a 1 MW FEL at X = 1 pm and Q = 5. Simple calculations based on (60)-(63) 
show that the mirror power density in a 20 m long cavity with z 0 = 0.2 is about 
1500 kW /cm 2 , which is quite large. However, if we reduce the Rayleigh length to 
z 0 = 0.02 and shorten the cavity length to 10 m, the mirror power density is reduced to 
about 600 kW/cm 2 . That is, we can reduce the power density on the mirrors by more 
than a factor of two and cut the cavity length in half by reducing the Rayleigh length by 
a factor of ten. 

Reducing the Rayleigh length is accomplished by making the mirrors larger and 
their radii of curvature smaller (recall z 0 = (2 SR - S 2 )' h l2L), allowing the light to cover 
a large mirror while focussed to a small waist at the center of the cavity. The stability 
cnteria for a symmetric optical cavity, 0 < (1 - S/f?)( 1 - S/ft) <1, reduces to z 0 2 2 0 
(Ref. 7). Therefore, as long as there is a Rayleigh length (i.e., z 0 > 0), the optical 
cavity will be stable. 

In section G of Chapter II, we discussed the physics of the FEL interaction in 
the Gaussian optical mode. We discussed how the optical field amplitude la G l is 
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greatest at t-x w , which causes the phase space separatrix height to increase as 
x : 0 i w then decrease as x : x w -> 1 . This effect, which is proportional to w”' /2 (x), 
complicates the phase space paths and will be exacerbated by a small Rayleigh 
length. The short Rayleigh length implies a small mode waist where the optical field 
amplitude will be much larger than at the ends of the undulator. We also saw how the 
phase of the Gaussian mode caused the separatrix to shift to the right as x : 0 -> 1 . 
This requires the initial phase velocity for maximum gain to be increased from v 0 = 2.6, 
the value found when using the plane wave approximation. For z 0 = 0.1, the phase 
shift for on-axis electrons is A<j> = 4>(x=1 ) - <)>(x=0) = 0.9n, and is smaller for off-axis 
electrons. With a small Rayleigh length, the wavefront curvature has to be larger in 
order that the light be focussed to a small waist and rapidly expand to cover the 
mirror. 

Operating the FEL with a short Rayleigh length is not the optimum design. The 
previous paragraph showed how the interaction of the electrons and the light changes 
with Rayleigh length. This concept has never been studied before now, precisely 
because the design is not optimum. In this chapter we will show the consequences of 
using a short Rayleigh length. Since we will analyze the case z 0 <1, we cannot 
expand (60) in 1/(z 0 ) as we did in Chapter II, but must use simulation instead. We will 
analyze the UVFEL of Chapter IV first, both the single pass gain G(x) and the gain 
evolution G(n) over several passes when in steady state. Next we present a MW- 
class FEL design based on the short Rayleigh length concept and simulations of 
same. 

B. UVFEL WITH SHORT RAYLEIGH LENGTH 

We first test the short Rayleigh length FEL concept by simulating the LPC 
UVFEL, which has parameters listed in Table 2. 
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1. Single Pass Gain 

We study the single pass gain G(x) in the UVFEL using transverse multimode 
simulations like those discussed in the previous chapter. In this case, however, since 
only one pass is simulated, cavity mirrors do not enter. The optical field is initialized 
as a Gaussian mode and can only evolve by mode distortion caused by the electron 
beam within the single pass. Contrast this with the multiple pass simulations 
presented in the previous chapter in which the optical field, entrained by the mirrors, 
evolves self-consistently over many passes. 

Plotted in Figure 22 is the single pass gain evolution G(-t) along the undulator 
for four different electron beam radii: (a) c e = 0.2, (b) c e = 0.45, (c) o e = 0.9, and (d) 
o e = 1.6. 




T 



Figure 22. Gain along the undulator G(x) for electron 
beams with radius (a) = 0.2, (b) c e = 0.45. (c) 

a e - 0.9, and (d) c e = 1.6. Rayleigh length is z 0 = 0.3, 
mode waist radius is w 0 = 0.5. 
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The Rayleigh length is z 0 = 0.3 giving a mode waist radius of w 0 = 0.5 (see (63)). The 
expected LPC UVFEL electron beam emittance e = 0.43 is again. Electron beams (a) 
and (b) have radii smaller than the mode waist radius (i.e., o e < w 0 ) so that they are 
enveloped by the light for the duration of their interaction. The radii of electron beams 
(c) and (d) are greater than the mode waist radius. Because of this, their final gain 
G(x=1) is reduced compared with beams (a) and (b). Also note that because of the 
large angular spread o e = 4.6 for electron beam (a), the final gain is smaller than that 
for beam (b). We saw in section D of the last chapter that electron beam (b) with 
c e = 0.45 and o e = 0.9 was nearly optimum for z 0 = 0.5 (optimum was o e = 0.4). 

The most encouraging aspect of the plots in Figure 22 is that although the final 
gain differs among the four electron beams, all show nearly exponential gain (see 
(58)) throughout the length of the undulator. This point is better illuminated by Figure 
23. In this figure we normalized the curves G(x) from Figure 22 to their respective final 
gains G(t=1) in order to compare the structure of the curves more closely. Figure 23 
shows convincingly that all four of the electron beams simulated resulted in similar 
gain evolution throughout the length of the undulator, even those ((c) and (d)) with 
o e > w 0 . We see that the gain evolution starts in all four beams after a small delay of 
t = 0.2, which is just the time it takes for the electron bunch to form (Ref. 3). 

In Figure 24 we have plotted G(t) for the same four electron beam radii, but 
with the Rayleigh length reduced to z 0 = 0.1. This means w 0 = 0.3 and only electron 
beam (a) with o e = 0.2 fits inside the optical mode throughout their interaction. For 
this reason, beam (a) has the highest final gain. Electron beam (b) has final gain 
reduced by nearly a factor of two over that in Figure 22 because it no longer fits inside 
the optical mode at the waist. Beams (c) and (d) also have final gain reduced, but 
only beam (d) has final gain less than 10%. This figure shows that in a single pass, 
where the electron beam cannot couple to higher order modes, it is more important to 
have the electron beam enveloped by the light than to have a small angular spread. 
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Figure 23. Normalized gain along the undulator 
G(t)/G(t=1) for electron beams with radii (a) o e = 0.2, (b) 
o e = 0.45, (c) o e = 0.9, and (d) o e = 1.6. Rayleigh 
length and mode waist radius same as in Figure 22. 



For Figure 25, the Rayleigh length was further shortened to z 0 = 0.02. In this 
case w 0 = 0.14 and none of the four electron beams fit inside the optical mode at the 
waist. The final gam for electron beams (a) and (b) is less than half their value at 
z 0 = 0.1 in Figure 24. However, there is still gain at this very short Rayleigh length, 
provided the electron beam is small enough. We see that c e = 0.9 or 1.6, curves (c) 
and (d), are clearly too large. The final gain G(t) > 20% for the two smallest electron 
beams. 
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Figure 24. Gain along the undulator G(x) for electron 
beams with radius (a) o e = 0.2, (b) a e = 0.45, (c) 
o e = 0.9, and (d) o e = 1.6. Rayleigh length is z 0 = 0.1, 
mode waist radius is w 0 = 0.3. 



In order to better compare the evolution of gain G(x) in Figure 25, the curves 
are normalized to their respective final values, as in Figure 23. The result is plotted in 
Figure 26. The effect of the short Rayleigh length is clear: single pass gain increases 
greatest where the fields la I are strongest, at the beam waist x = 0.5. The short 
Rayleigh length causes the fields at x = 0.5 to be much stronger than at any other 
position. This can be seen in electron beams (a), with o e = 0.2. and (d), with 
a e = 1.6. Electron beam (a) is focussed at x = 0.5; it was designed to take advantage 
of the field strength at the center by concentrating electron beam current there. 
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Electron beam (d) has a radius eight times that of beam (a). It is so large that gain is 
poor throughout the interaction at this Rayleigh length. What little final gain exists, 
however, was enhanced at x = 0.5 where optical fields are strongest. Indeed, we see 
in Figure 26 that there was net absorption for beam (d) until the field strength started 
to increase at x = 0.4. Also evident in the figure is that electron beams (b) and (c) 
have radii which allow gain to be distributed more uniformly. That is, they are larger 
than beam (a) so gain does not increase as much in the center of the undulator, but 
are smaller than beam (d) so that gain increases more at the ends of the undulator. 




x 



Figure 25. Gain along the undulator G(x) for electron 
beams with radius (a) o e = 0.2, (b) c e = 0.45, (c) 
Of, = 0.9, and (d) o e = 1.6. Rayleigh length is z 0 = 0.02. 
mode waist radius is w 0 = 0.14. 
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Figure 26. Normalized gain along the undulator 
G(t)/G(t = 1) for electron beams with radius (a) o e = 0.2, 
(b) c e = 0.45, (c) a e = 0.9, and (d) o e = 1.6. Rayleigh 
length and mode waist radius same as in Figure 25. 



To summarize, we have found that using a very short Rayleigh length degrades 
gam. However, gain is still significant for the LPC UVFEL at a Rayleigh length as 
short as z 0 = 0.02. For o e = 0.2, the final gain is G = 37%, while for a e = 0.45, the 
final gam is G = 23%. The gain is only senously degraded for electron beams with 
a e = 0.9 and o e = 1.6, which have radii two and four times greater than optimum, 
respectively. We also found that the most important region for the interaction is the 
center of the optical mode where the fields are strongest, even if the electron beam is 
much larger than the mode waist. This makes sense because it is the field amplitude 
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la I which drives bunching in the pendulum equation (36). 

When this research was started, however, it was thought that the interaction 
may "turn off" at the center of the undulator, meaning that section of undulator could 
be replaced by a drift space, making an optical klystron, which would improve gain. 
The optical klystron consists of two undulator sections with a drift space at the center 
(Ref. 3). The first section, called the modulator, starts the electron bunching process. 
The drift space lets the electron momenta continue bunching in the absense of the 
undulator field. The second undulator section, called the radiator, causes the bunched 
electrons to radiate coherently. Instead of a drift space, the same effect can be 
obtained in a shorter distance with a dispersive section, which is just a magnet that 
bends the electron beam out of and back into the beam path. Slower electrons bend 
through a smaller radius than the more energetic electrons. Since the faster electrons 
travel a greater distance through the dispersive section than the slower electrons, they 
are bunched upon exiting. It was found, however, that we do not need to resort to 
exotic klystron undulator designs to get the desired gain. 

2. Steady State Gain 

Having studied how the gain G(t) evolves within a single pass, we will now 
simulate several passes and follow the gain evolution G(n). Recall that G(n) follows 
the final gain G at x = 1 after each pass n. We use the transverse multimode, 
multiple pass simulations used in the previous chapter (see, for example, Figure 20) 
where the optical mode evolves self-consistently in the optical cavity over many 
passes. 

The electron beams (a), (b). and (c) from the previous figures were simulated. 
The simulations start from weak fields a 0 - 0.01 and evolve until the limit of weak 
fields la I = n is reached. At this point the final gain is in steady-state, by which we 
mean that single pass gain G(t=1) has nearly the same value over several passes 
and G(n) is nearly constant. It is this steady-state value of the gain evolution G(n) 
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that is plotted versus the Rayleigh length z 0 in Figure 27. 




Figure 27. Steady-state gain G versus Rayleigh length 
for electron beams with radius (a) o e = 0.2, (b) o e = 0.45, 
and (c) c e = 0.9. We find G > 60% at z 0 = 0.01 for 
beams (a) and (b). 



In Figure 27, we see that the nearly optimum electron beam, (b), has the 
highest steady-state gain at all Rayleigh lengths explored. This is in contrast with the 
single pass studies in which electron beam (a) had highest final gain for z 0 < 0.3. 
This electron beam is evidently more efficient in coupling to higher order modes than 
any of the other beams explored. Perhaps this is because of the better balance 
between it's beam radius o e = 0.45 and angular spread a 0 = 0.9 in their contribution to 
the electron phase velocity spread (i.e., it is better ’’matched'' (Ref. 3)). The maximum 
gam (or all three of the beams occurs at Rayleigh length z 0 = 0.1. The most important 
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feature to note in the figure shows the goal of this chapter: steady-state gain is 
substantial even at z 0 = 0.01. For electron beams (a) and (b) we have G(n)> 60%, 
while that for beam (c), twice the optimum radius, is still a respectable 12%. 

Our conclusion is that the FEL with a short Rayleigh length may be able to 
operate above threshold with respectable gain while spreading the optical mode at the 
mirrors. The trade-off is slightly lower gain for a much shorter optical cavity with 
reduced power loading at the mirrors. 

C. MW-CLASS FEL WITH SHORT RAYLEIGH LENGTH 

Encouraged by success in modelling the UVFEL with a short Rayleigh length, 
the concept was used to design a MW-class FEL. 



1. MW-Class FEL Design 

Our goal in designing a MW-class FEL is to derive a set of parameters which 
will provide sufficient gain within an optical cavity of length 10 m. This length was 
chosen because it fits within the beam of a modern naval cruiser or destroyer. 

FEL gain, discussed in section E of Chapter II, is dependent upon the 
dimensionless current density j (see (57) and (58)). Maximizing gain therefore 
requires that we maximize j. Recall that the definition of the current density in the 
linear undulator is 



y = 



BN (enK(J 0 (\) - MQ)L) 2 P 

Y*mc 2 



where % = K 2 /( 2(1 + K 2 )). Since N = L/Xq we can re-write (69) as 



j = 87T 2 - 



mc c 



K%JqW - Ji(S)) 2 l 3 



Y 3 



(69) 



(70) 



where we have separated the parts of j according to their source. Following the 
fundamental constants is the undulator contribution K 2 (JJ) 2 L 3 /X 0 where we have 
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shortened the Bessel function term in (70) to JJ. Since K ~ 1 in an FEL, which also 
fixes the Bessel terms, we are left with L and ^. 0 as designable undulator parameters 
(Ref. 3). The electron beam contribution is p/y 3 . The charge density p depends on the 
peak current t through 



P = 



Q 

necr?(l e lc) 



(71) 



where Q = 1 x (l e /c) is the charge per micropulse (not the cavity 
the electon beam radius, and l e is the length of the micropulse, 
energy y me 2 determines y. Therefore we can write 

. L 3 / 1 
' “ r| y 3 • 



quality factor), r e is 
The electron beam 
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Recall that the emission wavelength is given by 

X = V ( ~ ■ (73) 

2 r 

With K = 1 this becomes X = X 0 /y*. This means that if we fix the wavelength X, then 
^. 0 and y are not independent. Therefore, the remaining free parameters are 
L, /, r e , Xq, and X. 

Our philosophy in choosing a design was to use a short undulator and use 
peak current to maximize The reasons for this are twofold. First, we are using an 
optical cavity 10 m long with a short Rayleigh length. The undulator has to be short in 
order to fit around the optical mode at the center yet not "scrape” any of the light, 
which would damage the magnets. Second, a short undulator means the FEL will be 
less susceptible to electron beam quality. This is because the gain spectrum 
bandwidth, which is proportional to 1/(2 N), increases for decreasing undulator length 
L = NXq. This was discussed in section E of Chapter II and is important to the design, 
because as we demand more current from the accelerator we expect the beam quality 
to suffer. 
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